Elasticity and Response in Nearly Isostatic Periodic Lattices 
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The square and kagome lattices with nearest neighbor springs ol spring constant k are isostatic 
with a number ol zero-frequency modes that scale with their perimeter. We analytically study the 
approach to this isostatic limit as the spring constant k! for next- neare st-neighbor bonds vanishes. 
We identity a characteristic frequency uj* ~ \/W and length I* ~ \Jkjk' for both lattices. The shear 
modulus C44 = k' ol the square lattice vanishes with k' , but that for the kagome lattice does not. 

PACS numbers: 63.70.+h, 83.80.Fg, 64.70.kj, 63.20.D- 



An argument advanced by Maxwell [l[ states that rigid 
assemblies of loose particles must have at least Z lso = 2n / 
contacts per particle, where n/ is the number of rele- 
vant degrees of freedom per particle. This result has 
been applied to many systems liA including network 
glasses [H [H , rigidity percolation [jUsl , /3-cristobalite Q, 
granular media d E3 , protein folding [ill, and elasticity 
in networks of semi-flexible polymers [12J ■ 

Zero-temperature packings of particles can exhibit 
transitions with increasing particle density from an un- 
jammed, disordered state with no inter-particle contacts 
to a jammed, disordered state with an average of Z c 
contacts per particle. These transitions are discontin- 
uous because the coordination number Z c must sat- 
isfy Maxwell's inequality. For the special case of fric- 
tionlcss spheres, the coordination number at the jam- 
ming transition, Z c , appears to be exactly the mini- 
mal, or isostatic, value needed for mechanical stability, 
Ziso [IH EH EH ■ The coincidence of the jamming transi- 
tion with the threshold for mechanical stability gives rise 
to special properties at the transition: the coordination 
number Z jumps discontinuously but as the transition is 
approached from the high-density side, the shear modu- 
lus vanishes and length and time scales diverge as power 
laws [H QUE Eft In addition, the vibrational prop- 
erties just above the transition are very different from 
those of typical solids @, El, El El, El- 

The behavior of the isostatic jamming transition in- 
vites comparison to second-order structural transitions, 
in which the vanishing of an elastic constant in an iso- 
static lattice signals the instability towards a shape dis- 
tortion. Jamming involves a transition between two dis- 
ordered states, while critical structural transitions in- 
volve a crystalline state. This raises the question of 
which aspects of the jamming transition apply to all sys- 
tems, periodic or disordered, near isostatic transitions 
and which do not. To answer these questions, it is impor- 
tant to study models for which results can be rigourously 
established via analytic calculations. In this paper, we 
test the robustness of the connection between isostatic- 
ity, power-law elastic moduli and vibrational properties 
by an exact analytical exploration of the approach to 
the threshold of mechanical stability in fully periodic, 
nearly isostatic systems. Specifically, we study two lat- 



FIG. 1: (Color Online) (a) Square and (b) kagome lattices 
with NN springs of spring-constant k and NNN springs of 
spring constant k' . White circles in (a) and white triangles in 
(b) show a zero-energy distortion. 



tices, the 2D square and kagome lattices, shown in Fig.[TJ 
with nearest-neighbor (NN) harmonic springs of spring 
constant k and next-nearest-neighbor (NNN) harmonic 
springs with spring constant k' . The isostatic structural 
transition is then approached continuously as k' — > 
since both systems are isostatic there with Z c = 4. To 
describe the phase that results when k' < 0, it is neces- 
sary to add a nonlinear term gx 4 to the energy [20| . 

Our principal results are that in both the square and 
kagome lattices, there is a characteristic frequency to* ~ 
y/k' and a length I* ~ \TkJk' that follow directly form 
the form of the wavenumbcr- and frequency-dependent 
response function but that can also be obtained from 
cutting arguments [5. ITtt However, the shear modulus 
exponent is different in the three systems, showing that 
isostaticity does not confer universality on all power-law 
properties near the transition. 

It is useful first to understand the two states that strad- 
dle the structural transition. Above the transition, where 
k' > 0, the square/kagome lattice is stable. Below the 
transition, where k' < (and g > 0), the A*o zero modes 
at k' = develop positive or negative amplitudes, and 
there are ~ a^ , where a > 1, distinct ground states. 
In this paper, we will restrict our attention to the case 
k' > 0, and we will use the harmonic approximation. 

Now consider the square lattice shown in Fig.QJa) with 



NNN nonlinear springs. This lattice has TV = N x N y 
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FIG. 2: (a) ^(q) when fc' = showing the line of zero energy 
at q x = and one- dimensional dispersion as a function of q x 
at fixed q y . (b) Density plot of the low-energy phonon mode 
when k' = 0.02fc showing cos4# modulation at small q and 
one- dimensional isostatic behavior at large q. 



bonds. Thus when k' = 0, this system is isostatic with 
a number of zero-frequency modes, No = 2N — N nn — 
N x +Ny, equal to half the perimeter of the system. When 
k! > 0, these modes become quasi-isostatic with nonzero 
frequencies that vanish as k' — » 0. The components of 
the dynamical matrix are easily calculated |21j 

D xx (q) = D yy (q y , q x ) = 4fcsin 2 (g x /2) + 4/c' sin 2 (g s/ /2) 

+ 4fc'sin 2 ( (?x /2) - 8k' sin 2 (g a /2) sm 2 (q y /2) 

D xy(<l) = D yx(q) = 2k' sin(fo) sinfaj,), (1) 

where we set the lattice constant a equal to 1. In the 
continuum limit ijCl, Z?y (q) obtains the form dictated 
by elasticity theory with elastic constants C\\ = C xxxx = 

Cyyyy, C\2 = C XX yy and C 44 = C x y x y'. [2ll | 



D xx (q) = C n ql + C 44 g 2 ; D yy (q) = Cxiqfr + C 44 q 2 
D x V (q) = (C12 + C 4 4)q x q v , (2) 

and we can relate k and fc' to the elastic constants: 
C\\ = k + k', C 44 = C\2 = k' . Thus, the shear modulus 
vanishes as k' — > 0. When fc' = 0, Dy(q) breaks up into 
two independent one-dimensional compressional phonon 
systems: it is diagonal with D xx (q) = 4fc sin 2 (q x /2) inde- 
pendent of q y and D yy (q) = 4fc' sin 2 (g y /2) independent 
of q x . Thus at q x = 0, -D Ka; (q) vanishes for all points 
along the line — n < q y < tt. Note that at a standard 
structural phase transition, components of Dy (q) vanish 
at one or possibly a discrete set of points in the Bril- 
louin Zone (BZ) when an elastic modulus vanishes. In 
contrast, for the periodic isostatic system, components 
of Dij (q) vanish along lines. 

The one-dimensional nature of Z?y(q) gives rise 
to compressional phonons with frequencies to x _ y (q) = 
2 yk'l shxq Xt y/2\ [Fig. [2{a)] and a one-dimensional den- 
sity of states p(u))(2/tt) / \/4k — u> 



2 with a nonzero value 



1 at oj = as shown in Fig. 03a). 

> 0, modes exhibit a cos 40 modulation 



po = (nVk) 
When k' 

at low frequency and one-dimensional isostatic behav- 
ior at larger q [Fig. Efb)]. When < k' < k, 
Dij(q) is well approximated as a diagonal matrix with 
D xx (q) = kq 2 + 4k' sm 2 (q y /2) with associated eigenfrc- 



quency uj x (q) ~ yJ ' D xx (q). These expressions immedi- 
ately define a characteristic frequency 

uj* = 2-Jk 1 (3) 

as the frequency at the point M = (0, ir) on the BZ edge. 
The first term in -Da^q) represents the long- wavelength 
anomalous isostatic (iso) modes that are present when 
k' = 0, whereas the second represents the effects of NNN 
coupling. When q x = 0, the only length scale in the prob- 
lem is the unit lattice spacing, and no divergent length 
scale can be extracted from D xx (0,q y ) as it can be in 
the case of the structural phase transition. When the 
first term is large compared to the second, -D^q) re- 
duces to its form for the isostatic k' = limit, and we 
can extract a length by comparing these two terms. The 
shortest length we can extract comes from comparing kq 2 
to -D^q) at point M on the zone edge, i.e. 



T = {1/2)^1^ = l/q* 



(4) 



If q y < 7T, the isostatic limit is reached when q x > q*. 
A similar analysis applies to D yy (q) when q y > q* . If a 
square of length I is cut from the bulk, the wavenumbers 
of its excitations will be greater than tt/1, and for ql* > 1, 
all modes within the box will be effectively isostatic ones. 
This construction is equivalent to the cutting argument 
of Wyart et al. [3,[22|. 

Equation (j4|) is identical to the length at which the fre- 
quency of the compressional mode 0J x (q x , 0) = \fkjV* ~ 
\JC\\jl* becomes equal tow*. A meaningful length from 
the transverse mode Lu x (0,q y ) cannot be extracted in a 
similar fashion. The full phonon spectrum [Fig. [3ja)] ex- 
hibits acoustic phonons identical to those of a standard 
square lattice at q <C 1 and a saddle point at the point M. 
Thus, the low-frequency density of states is Debye-like: 
p(uj) = (uj I \2tt)) I \J kk' with a denominator that, because 
of the anisotropy of the square lattice, is proportional to 
the geometric mean of longitudinal and transverse sound 
velocities rather than to a single velocity. In addition 
p(uj) exhibits a logarithmic van Hove singularity at oj* 
and approaches the one-dimensional limit (\/ir)/yfk at 
w* « w « 2\fk. The frequency uj* [Eq. is recovered 
by equating the low-frequency Debye form at uj* to the 
high-frequency isostatic form of the density of states. 

The kagome lattice can be viewed as an array of one- 
dimensional staggered linear rows, parallel to the x-axis, 
of pairs of opposing triangles as shown in Fig. [TJ Iden- 
tical arrays with rows parallel to (cos 27r/3, ± sin27r/3) 
can be identified. As in the square lattice, each site has 
2c? = 4 nearest neighbors in the bulk. A counting proce- 
dure similar to that used for the square lattice yields Nq 
proportional to the s/~N for lattice of N sites. 

The isostatic (iso) modes correspond to identical ro- 
tations about their top vertices of all of the "up" point- 
ing triangles in any row in the horizontal (or symmetry- 
equivalent) grid. These rotations require counter rota- 
tions of connected "down" pointing triangles as shown 
in Fig. [lib). There are three sites per unit cell (which 
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FIG. 3: (Color Online) Density of states as a function of 
r\ = u)/uj* for (a) the square and (b) kagome lattices, showing 
their constant value when k' = (dashed green line) and van 
Hove singularities, from the saddle at S and minima at M 
and r, whose frequencies all scale as ui* oc vF when k! > 0. 
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FIG. 4: (Color Online) (a) Phonon dispersion along symmetry 
directions. The dotted lines are for k 1 = and the solid lines 
are for k' = 0.02. The isostatic and quasi-isostatic branches 
are in red. (b) shows isostatic and shear modes along YM 
and indicates characteristic frequencies and wavenumbers. 



we take to be those in the "up" triangles) in the kagome 
lattice and six phonon branches [Fig. 0Ja)]. Three of 
these are high-energy optical branches, two are acous- 
tic and one is isostatic. The zero modes of the latter 
show up as three lines of zero frequency along q x = 
[r = (0,0) to M = (0,G /2), where G = 47r/V3, in 
the BZ] and the two symmetry-related lines as shown in 
Fig. [5]Ja). Away from q x = 0, the isostatic mode fre- 
quency is w/(q) = c x q x , where c x = V~3k/4 for small q x . 
This behavior is identical to that of the square lattice. 
The resulting density of states decreases linearly with oj 
from a nonzero value po — 2>Gq/ {Air 2 c x ) = 8/(n\/k) at 
u> = 0. The total low-frequency density of states from 
the two acoustic modes and the isostatic mode is inde- 
pendent of oj and equal to po at small lu [Fig. |3fb)]. 

When springs of spring constant k' are added to the 
NNN bonds shown in Fig.fljb), the quasi-isostatic mode 
along q x = [Fig. IDJb)] has nonzero frequency of order 
s/W for all q y > and gives rise to various lengths of order 
/* ~ yjk/k'. At q x = and low values of q y , this mode 
hybridizes with the transverse phonon mode to produce 
a gapped translation-rotation mode with frequency 



to 



'6k' 



(5) 



at q = 0. At small q, there are isotropic longitudi- 
nal and transverse sound modes with respective veloc- 



FIG. 5: (a) Plot of iso mode at k' = 0. (b) density plot of iso 
mode at k' = 0.02fc. Inset: circular symmetry at origin. 



ities cl = \/3k/4 and ct = v / fc/4. The lowest fre- 
quency mode after hybridization at q x = is a transverse 
phonon near q y = and predominantly an isostatic ro- 
tation mode at q y > q* H , where q* H = A^j2>k' jk = l/l* H 
can be termed a hybridization wave number. This mode 
reaches a maximum frequency lo* s = uo* /\/2 at a sad- 
dle point at q y = Q s = 4(3fc'/2fc) 1 / 4 - (P)- 1 / 2 and a 
local minimum with frequency uj\j = u>* at the zone- 
edge point M. At low frequency, the DOS is Dcbyc-likc: 
p(uj) = (uj/2tt)(c~ l 2 + c T 2 ) = 32uj/(2nk). The points M, 
r and S = (0, Qs) give rise to van Hove singularities 
in the DOS. The minimum point A/ produces a jump 
ApM = 8V2/n\/k > po at the saddle S produces a 
logarithmic singularity at wj, and the minimum at T a 
jump Apr = (16/7r)^/3fc'/2fc <C po, which is just visible 
in Fig. [!». 

Lengths that scale as \fk~jkf can be introduced in much 
the same way as in the square lattice. The square of the 
low-frequency isostatic-shear mode at q x = increases 
as c x q x for nonzero shear, and lengths 2< 
and I 



(1/4)VW 

M — (v / 6/8)yfc/fc' follow from comparing c x q x to 



and 



M ' 



respectively. These two lengths are longer 



than the hybridization length, l* H . Thus, it is only at 
length scales less than l* H that isostatic modes are re- 
trieved completely, and we should take 



r = r H = (i/4)7fc7(3fco. 



The less divergent length £s = Q c 



(6) 

I* determines 

the position of the saddle S and does not describe the 
same physics as the other lengths. 

The long-wavelength, low-frequency properties of the 
kagome lattice are best understood by considering the 
effective low-energy dynamical matrix Dij obtained by 
integrating out the high-energy modes. This matrix 
is conveniently represented in a basis consisting of the 
longitudinal (L) and transverse (T) phonons and the 
mixed rotation-translation gapped mode (77) at q: Dlt = 
Dtl = 0, and components are 

3kq 2 /16; D TT = kq 2 /16; D m = 6k' + kq 2 /16 



Dll 

D Ln = D* L = kq 2 cos(30)/16 (7) 

D Tr) = D* T = -i(V3/2)k'q + kq 2 sm(W)/l6, 

where 6 is the angle that q makes with the x-axis. When 
k' = 0, the transverse phonon and rotation-translation 
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TABLE I: Frequencies, lengths, and elastic moduli in square 
and kagome lattices and the marginally-jammed (MJ) state 
of disordered sphere packings, u* , {* and B scale the same 
way in all lattices if we take k! ~ (Az) 2 , but G does not. 



Quantity 


Square 


Kagome 


MJ 


r 

UJ* 

G 


(fc') i/a 
k' 


k 


(Az)- 1 
(Az) 1 
(Az) 1 


B 


k 


k 


(Az)° 



mode mix at q x = (6 = n/2) to produce the zero 
isostatic rotation mode and a sound mode with to — 
yjk/$q y . In addition, there is the longitudinal mode 
with u> = (y3k/A)q y . When k' ^ 0, the off-diagonal 
terms can be ignored to lowest order in q, leading to 
longitudinal and transverse sound modes with respective 
velocities cl = yB and ct = VG, arising from bulk and 
shear moduli B = 3fc/16 and G = fc/16. Note that the 
shear modulus is proportional to k; this is in contrast to 
the square lattice where G oc k' . To understand how this 
comes about, it is instructive to look at the TT compo- 
nent of the phonon susceptibility, 



XttM) = D T ^(q,9) = 



16 6 + (g/*) 2 (3 



COS' 



36») 



kq 2 6 + 2(ql*) 2 cos 2 3(9 



(8) 

When ql* <C 1, this reduces to the required isotropic form 
l/(Gq 2 ). When ql* ^> 1, xtt is anisotropic with value 
l/(Gg 2 ) for cos 36* = 1 and l/(6fc') for cos 36* = 0. 

We can now compare the properties of nearly iso- 
static lattices with those of disordered sphere pack- 
ings. Marginally-jammed packings with volume frac- 



tion (f> = (fi c + A(p, just above the volume fraction <j) c 
at the jamming threshold, have an average number of 
contacts per sphere z = z c + Az, where z c = 2d, 
and Az ~ (A0) 1 / 2 [l3|, They are macroscopically 
isotropic with bulk and shear moduli that for harmonic 
inter-particle potentials scale respectively as B ~ (Az)° 
and G ~ (Az) 1 [H, [TB|. The longitudinal and trans- 
verse sound velocities then scale as cl ~ (Az)° and ct ~ 
(Az) 1 / 2 . The density of states has a plateau above a fre- 
quency uj* [H, EH- A length scale I* ~ 1/Az can be ex- 
tracted by equating c^/l* to uj* or from the response 
to a point perturbation [23|. This scaling emerges from 
the cutting arguments of Refs. H, [l?]- A second length 
I*, (Az)" 1 / 2 can be defined from cr/l^r = oj* [H] and 
is important for energy transport [13]. 

The most robust features are the scalings of uj* ~ y/k* 
and I* ~ in both the square and kagome lat- 

tices. Effective medium calculations [24| for these lattices 
with NNN bonds added with probability P ~ Az yield 
k' ~ (Az) 2 ; this correspondence yields the observed scal- 
ings of uj* and I* for marginally-jammed systems. These 
scalings are likely to be robust for all nearly-isostatic 
systems. However, the power-law scaling of the shear 
modulus is not universal (Table [T|. The shear moduli of 
marginally-jammed systems and the square lattice both 
vanish with Az (with different powers), but G of the 
kagome lattice is proportional to k and does not van- 
ish. Thus, different isostatic transitions can fall into dif- 
ferent universality classes. This conclusion is consistent 
with recent numerical results for disordered systems near 
isostaticity, which show different scalings of bulk mod- 
uli m. 
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